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Abstract 

The relativistic Hartree approach describing the bound states of both nucleons 
and anti-nucleons in finite nuclei has been extended to include tensor couplings for 
the u- and p-meson. After readjusting the parameters of the model to the properties 
of spherical nuclei, the effect of tensor-coupling terms rises the spin-orbit force by 
a factor of 2, while a large effective nucleon mass m*/M^ ~ 0.8 sustains. The 
overall nucleon spectra of shell-model states are improved evidently. The predicted 
anti-nucleon spectra in the vacuum are deepened about 20 - 30 MeV. 
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I. INTRODUCTION 


One of the main characters distinguishing relativistic approaches from nonrelativistic 
approaches is that the former one has a vacuum. It is quite interesting to study the 
structure of quantum vacuum in a many-body system, e.g., in a hnite nucleus where the 
Fermi sea is hlled with the valence nucleons while the Dirac sea is full of the nucleon-anti¬ 
nucleon pairs. In the relativistic treatment of nuclear phenomena [1], the Dirac equation 
is used to describe the behavior of nucleons in nuclei. The effects of the nuclear medium 
on nucleons are taken into account through introducing strong Lorentz scalar (S') and 
time-component Lorentz vector {V) potential. In the language of meson-exchange theory 
the scalar potential can be attributed to the exchange of sigma meson and the vector 
potential to the exchange of omega and rho mesons as well as the electromagnetic force. 
Since the Dirac equation describes the nucleon and the anti-nucleon simultaneously, the 
effects of mean helds act on both of them. Consequently, not only the valence nucleons are 
bounded in the shell-model like states, but also there exist bound states for anti-nucleons 
emerging from the lower continuum. The observation of anti-nucleon bound states is a 
verihcation for the application of the Dirac phenomenology to a relativistic many-body 
system [2]. It constitutes a basis for the widely used relativistic mean-held (RMF) theory 
[1, 3, 4, 5, 6, 7] and the relativistic Hartree approach (RHA) [8, 9, 10, 11, 12]. Since 
the bound states of nucleons are subject to the cancellation of two potentials S + V {V 
is positive, S is negative) while the bound states of anti-nucleons, due to the G-parity, 
are sensitive to the sum of them S — V, consistent studies of both the nucleon and the 
anti-nucleon bound states can determine the individual S and V. In addition, the exact 
knowledge of potential depth for anti-nucleons in the medium is a prerequisite for the 
study of anti-matter and anti-nuclei in relativistic heavy-ion collisions [13, 14]. 

The shell-model states have been theoretically and experimentally well established [15] 
while no information for the bound states of anti-nucleons in the Dirac sea are available. 
This is the aim of our work. In Ref. [16] we have developed a relativistic Hartree ap¬ 
proach which describes the bound states of nucleons and anti-nucleons consistently. The 
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contributions of the Dirac sea to the source terms of meson-held equations are considered 
up to the one-nucleon loop and one-meson loop and evaluated by means of the derivative 
expansion technique [17]. The parameters of the model are adjusted by htting to the 
properties of spherical nuclei. The major outcome of the RHA model is that a rather 
large effective nucleon mass m*/M n ^ 0.8 is obtained compared to the value of 0.6 in 
the relativistic mean-held calculations where the no-sea approximation is adopted . This 
is caused by the ehects of the vacuum contributions which decrease the magnitude of 
the scalar potential S substantially. Correspondingly, the vector potential V is also sup¬ 
pressed since the quantity R -|- S' is controlled by the saturation properties. As pointed 
out above, the anti-nucleon bound states are mainly determined by the sum of the scalar 
and vector potentials S — V. The smaller values of S and V obtained in the RHA cal¬ 
culations lead to a weaker bound on the single-particle energies of anti-nucleons, which 
turn out to be only half of that computed in the RMF model. On the other hand, the 
spin-orbit potential of nucleons is related to d{S — V)/dr. In the RHA calculations the 
spin-orbit splitting of the shell-model states is roughly 1/3 of that calculated in the RMF 
approach and indicated by the empirical data, although the general trend of the energy 
spectra coincides with each other. Because our goal is to develop a model to predict the 
bound states of anti-nucleons in the vacuum with the model parameters constrained by 
the nuclear bulk properties, in order to get reliable results for the anti-nucleon spectra 
one should first describe the nucleon spectra as good as possible. 

Theoretically one can incorporate tensor couplings for the u- and p-meson, which 
mainly contribute to the spin-orbit force. The model now becomes non-renormalizable. 
However, from the point of view of modern effective held theory [18] the argument of 
renormalizability is not a severe restriction to theories. Since we will employ an effective 
Lagrangian of mesonic degrees of freedom and nucleons, just as the Skyrme force is an 
effective Lagrangian for nonrelativistic calculations [19], the mesons here are effective 
mesons. In an effective held theory one normally evaluates the lowest-order diagrams and 
makes regularization whenever a divergence appears. The parameters of the Lagrangian 
are adjusted to ht certain experimental data. The validity of the whole approach is 
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justified by successful explanations and predictions of observables. In other words, within 
the framework of an effective held theory one mainly concerns the balance between the 
predictive ability of the theory and the complexity of the theory. 

In the present work we will investigate the effects of tensor couplings in the relativistic 
Hartree approach for hnite nuclei. In this extended version of the model the parameters 
will be rearranged in a least-sqnare ht to the properties of spherical nnclei. The model 
is then applied to stndy the bonnd states of nucleons and anti-nncleons. The paper is 
organized as follows: In Sect. II we introdnce the effective Lagrangian and review the 
RHA model. In Sect. Ill we present the nnmerical results and discussions. A snmmary 
and outlook are dually given in Sect. IV. 

II. RELATIVISTIC HARTREE APPROACH 

The Lagrangian density of nucleons interacting through the exchange of mesons can 
be expressed as [1] 

C = C-p + Cl- ( 1 ) 

Here is the Lagrangian density for free nucleon, mesons and photon 

+ ^di^ad^a -U(a)- 

+ - jB.,, ■ R'“ + ■ R" -(2) 

and U((t) is the self-interaction part of the scalar held [20] 

U (a) = ^rnla"^ + (3) 

In the above expressions ip is the Dirac spinor of the nucleon; a, R^ and represent 
the scalar meson, vector meson, isovector-vector meson held and the electromagnetic held, 
respectively. Here the held tensors for the omega, rho and photon are given in terms of 
their potentials by 



(4) 


(5) 

- d^A^. 

(6) 
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Cl is the interaction Lagrangian density 

f _ 1 - 

Cl = - -gp'ipl^.r ■ 

4i\a_/v ^ 

- ^e^(l + TQ)-i^^A^. (7) 

Here = | [ 7 ^, 7 ^], r is the isospin operator of the nucleon and tq is its third component, 
go-, g^, gp and jA.'n = 1/137 are the coupling strengths for the a-, a;-, p-meson and for 
the photon, respectively, and fp are the tensor-coupling strengths of vector mesons. 
Mtv is the free nucleon mass and nia, rriuj, rup are the masses of the a-, a;-, and p-meson. 
In hnite nuclei the Dirac equation is written as 
d ~ f 

■ V+ /3(Mjv-gac^(x))+g<^o;o(x) - ^^*7 ■ (Vo;o(x)) 

+ ^gpT-o^o,o(x) - ^^fro 7 ■ (Vi?o,o(x)) + ^e(l+ro)Ho(x) ^/(x,t).(8) 

The field operator can be expanded according to nucleons and anti-nucleons and reads as 


^(x, t) = J2 [bai’a{^)e . (9) 

a 

Here the label a denotes the full set of single-particle quantum numbers. t/’a(x) are the 
wave functions of nucleons and t/’a(x) are those of anti-nucleons; and are their 
positive energies, respectively, ba and d'^ are the annihilation and creation operators 
of nucleons and anti-nucleons that satisfy the standard anticommutation relations. We 
assume that the meson helds depend only on the radius and discuss the problem in 
spherically symmetric nuclei. In this case, the usual angular momentum and parity are 
good quantum numbers. As described in Refs. [ 21 , 22 ], the eigenfunctions of nucleons 
are the well-known spherical spinors 


^a(x) = 




We make the ansatz for the wave functions of anti-nucleons [16] 



( 10 ) 


( 11 ) 
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Here Vtjim are the spherical spinors dehned as 


^ I m'msW^ Yim'Xlras^ 


Yim' are the spherical harmonics and x^-ms eigenfunctions of the spin operators. 

Ga, Fa and Fa, Ga are the remaining real radial wave functions of nucleons and anti¬ 
nucleons for upper and lower components, respectively. 

Inserting Eq. (9) into Eq. (8) one immediately obtains two relativistic wave equations 
for the and '^“(x). Applying the concrete expressions of the wave functions given in 

(10) and (11), we arrive at the coupled equations for the radial wave functions of nucleons 


EaGa{r) = --^ + — -7^idrMr))-j^roa{drRo,o{r)) Fair) 

dr r 2 Mn 4:Mn 

r 1 

Mjv - gaO-(r) gujUJoir) + -gpToaRofiir) 

-|--e (1 -|-Toq) Ao(r) ^^(r), (13) 

EaFair) = ^ idrUJo{r)) - T^roa (drRoflir)) Ga(r) 

dr r 2M]sf 

+ -Mm + gaO-(r) g^ujQ^r) + ]^gpToaRofiir) 

+ ^e(l-hroa) Ao(r) Fair) (14) 


and anti-nucleons 


-EaFair) — idr^oir)) + -^^ToaidrRofiir)) Gair) 

dr r 2 Mm AMm 

+ Mm - gacrir) + gu^uJoir) + ^gpfoaRoM 

-F^e(l- toq) Ao(r) Fair), (15) 

-EaGair) = -^ + — + ^^idr^oir)) + ^^fQaidrRofiir)) Fair) 

dr r 2 Mm AMm 

+ -Mm + ga(^ir) + gujUioir) + ^gpfoaRofii'f') 

-F^e(l- toq) Ao(r) Gair), (16) 
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where 


[/ forj=/-i 

and toq is the isospin factor of anti-nncleons, foa = —tqq. In the nnmerical solntion of the 
relativistic wave equations one eliminates the small components to obtain the Schrodinger- 
equivalent equations. For the nucleon we eliminate the lower component while for the 
anti-nucleon the upper component. By dehning the Schrodinger-equivalent effective mass 
and potentials of the nucleon 

Meff = Ea + Mm - gaO-(r) - gujUJo{r) - ]^gpToaRo,o{r) - + roa) Ao(r), (18) 

Ueff = Mm- gaCr{r) + gu;UJo{r) + ^gpToaRofi{r) + ]^e{l + Toq) (19) 

^(^) = y - {drUjQ{r)) - {drRofi{r)) (20) 

and the anti-nucleon 


Meff = Ee, + Mm- gaCr{r) + g^Uoir) - ^gpro„i?o,o(^^) + ^e{l + ro„) Ao{r), (21) 
Ueff = Mm- gaO-(r) - ga;Cc;o(r) -F ]^gpToe,R^fi{r) - ]^e{l + Toa) Ao{r), (22) 

W(r) = - + (drRiAr)) . (23) 

r 2Mm AMm 

we arrive at the Schrodinger-equivalent equations for the upper component of the nucleon’s 
wave function 


EaGa{r) = 




M. 


-1 

eff 




Gai'f') + UeffGa{r) 


and the lower component of the anti-nucleon’s wave function 


EaGa{r) = 


1 + »•(,') 


M. 


-1 

eff 




Ga (r) + UeffGa {r ). 


The small components can be obtained through the following relations 

d 


F„(r) = 


dr 


+ W{r) 


Eair) = -M^ff 


d 


'^ + IT(r) 

dr 


G«(r), 
Ga{r). 


(24) 


(25) 


(26) 

(27) 
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In the above we have changed toq —>■ —Toa, i.e., now the anti-nucleon has the same 
isospin factor as the corresponding nucleon. From Eqs. (24) and (25) one hnds that the 
Schrodinger equation of the anti-nucleon has the same form as that of the nucleon. The 
only difference relies on the dehnition of the effective mass and potentials, that is, the 
vector helds change their signs. The so-called G-parity comes out automatically. 

It’s now a suitable place to discuss the orthonormalization condition of the wave func¬ 
tions. In Eq. (9) the annihilation and creation operators of nucleons and anti-nucleons 
satisfy the usual anticommutation relations. From the equal-time anticommutation rela¬ 
tion of the Dirac held operator one can derive the matrix equation for the normalization 
of the wave functions 

^ [^a(x,t)^/’+(y,t)+<(x,tX+(y,t)] =5(x-y). (28) 

a 

In static case one can eliminate the t index. The following orthogonal conditions can be 
obtained from the Dirac wave equations [18] 

J = 0 ii a (3, (29) 

J = 0 if a 7 ^/3, (30) 

y d^a;^/>+(x)^/>^(x) = 0. (31) 

Multiplying Eq. (28) on the right with t/’a(y) or t)e found that these wave 

functions must satisfy the normalization conditions 

J A^^(y)^a(y) = J {yWaiy) = (32) 

This leads to the orthonormalization conditions for the radial wave functions of nucleons 
and anti-nucleons 

fOO 

/ dr [Gc,{r)G/3{r) + Fo(r)F^(r)] = 6^,0, (33) 

J 0 

r — — — — n 

/ dr Gc,{r)Gp{r) + Fo(r)F^(r) = 5 ^ 0 , (34) 

Jo L J 

respectively. The single-particle energies of the nucleon and the anti-nucleon can be 
evaluated as 

coo \ d 1 \ d 

= / dr{G^{r) -- + W{r) F„(r)+F„(r) - + W{r) G^{r) 

Jo I dr J [dr 
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( 35 ) 


+G„{r)U,f,G„{r) - F„{r) [M,„ - B„] F„(r)}, 

Eo = / cir{-Fo(r) 

Jo 

+Ga{r)UeffGa{ 

which are obtained through the iteration procedure. 

The main ingredients in Eqs. (18) - (23) are the meson helds, which are determined 


by the Laplace equations 

(v" - K] 

1 1 

1 (rir) = -gaPsir) + ^ba^ir) + —C(T^(r), 

(37) 

1 

- ml] 

)^o(r)= guppoir) 

(38) 

1 

- ml) 

Fo,oir) = ^gpPofiir) 

(39) 

VMo(r) = 

—^PPrfiir), 

(40) 

and 


cP 2d 

V =-^- 

(41) 


in the case of spherical nuclei. The source terms of the meson-held equations are various 
densities which, in principle, should contain the contributions both from the Fermi sea 
and the Dirac sea. Under the mean-held approximation, those densities are formally 
computed as the expectation values of various bilinear products of held operators in the 
ground state of the many-body system. For example, a direct calculation of the baryon 
density gives 

po = {i’o I i’*i> I i’r>) 

= E (42) 

a=l Q:=vac 

where the sum on the second term of the last equality runs over all anti-nucleon spectra 
in the vacuum, and therefore causes divergence. A proper regularization scheme is appar¬ 
ently needed in order to render it to a hnite value. Unfortunately, at the situation of hnite 
nuclei it is currently untractable. If one simply cuts oh it, the integration of the term 
with respect to the space doesn’t vanish. This violets the baryon number conservation. In 


d - , , 



d - , , 


G„(r) 

-G„(r) 



Fair) 


Fair) 


M, 


eff 




Fair)}, 


(36) 
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an alternative point of view, the mean fields are taken as a starting point for calculating 
corrections within the framework of Quantum Hadrodynamics [23]. That is, one drops 
the second term and includes quantum corrections by means of Feynman diagrams and 
path-integral methods. Within this scheme, the contributions of the valence nucleons to 
the densities are computed by adding up wave functions while the contributions of the 
Dirac sea are taken into account in loop expansions. At the one-loop level the effective 
action of the system can be written as [24] 

+ Ct'^ + r^aience + ^^Tr lu - ^hTr In . (43) 

Here the CT are the counterterms. Fvaience is the contribution from the valence nucleons, 
which for time-independent background helds is just minus the energy of the valence 
nucleons. The last two terms in Eq. (43) represent the contributions of the Dirac sea 
stemming from the one-meson loop and one-nucleon loop, respectively. By means of the 
derivative expansion technique [17] they can be expressed as 

r«(cT) = lhTrln(iD-^) 

= J + , (44) 

rW(^) = -ihTrln{iG-^) 



It can be verihed that the tensor-coupling terms contribute to the higher-order terms in 
the derivative expansion and have therefore been neglected. vj^\a) and Vp^\a) are the 
effective potentials from the one-meson loop and one-nucleon loop, in which the field is a 


constant, o'(a;) = do, the same situation as in nuclear matter. These two terms contain 
divergent part and should be regularized. Through adding suitable counterterms, {a) 
and Vp\a) can be calculated in nuclear matter which turn out to be [25, 26, 1] 



m* / ba ca'^ \ ^ / ba ca‘^ \ 

(87r)2 y ml 2ml J ^ ml 2ml J 


/ ba ca^ \ 
[ml ^ 2ml) 
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3/6(7 C(7^ \ ^ 1 / ba\^ / ba 3c(7^\ 1 / ba\^ , , 

2 \ml 2mlJ 3 \ml) \ml 2ml / "*” 12 \ml J 

Vj^\a) = (Mjv - ga(T)^ In (^1 - + M^g^a - ^M^gla^ 

+yM7vg^a^ - . (47) 


The functional coefficients before various derivative terms can be determined in the deriva¬ 
tive expansion technique and read as [9, 12] 


z<‘V) 

diV) 

dl(ff) 

dlV) 


1 ( 6 - 1 - caY 

12 lQ7r‘^{ml + ba + Ica"^) ’ 






-^In 

37r^ 







(48) 

(49) 

(50) 

(51) 


Inserting Eqs. (44) - (51) into Eq. (43) and minimizing the effective action with respect 
to the corresponding helds, one reveals the meson-field equations as given before, in which 
Tvaience coustitutcs the densities originated from the valence nucleons while r(^)((7) and 
r*^^)('/’) compose the densities stemming from the Dirac sea. At the end we obtain concrete 
expressions of various densities contained in Eqs. (37) - (40) 


Ps{r)=pf{r) + pr{r), 

(52) 

po(r) = prw+prw, 

(53) 

Po(0 = Po’^“'(r), 

(54) 

po,o{r) = Plpir), 

(55) 

ploir)=P^o’r\r), 

(56) 

pPr,oir) = Pp“'o(^) + Pp“o(^)i 

(57) 


and 



1 

47rr^ 


Y,w^{2j^ + 1) [Gl{r) 


a=l 



(58) 
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pT\r) = 


47rr2 


Wa{2ja + 1) G 


I] w^a(2ja + 1)2 [(9^G'c,(r)F«(r)], 

Q:=l 

Poro(^) = H M^a(2ja + l)roa [G'^(r) + F2(r 

Q:=l 


T.val( \ 

Po.o = 


-—r ^ Wci2ja + l)roa2 [i9^Ga(r)F„(r)], 
4:7rG 


pSoW = ^[p^W + pSW], 


pr(o = 


K’w ++ y f^r 

(90- L J 2g<^ dc 


Z(«(ff)l - 


pr(^) = -S 


g.5aL'« 2g<, [da \ [dr ^ 

, 1 v(l)^ ^ , 2 C? ^ ^ ^ , ga r ^ ^ 

_ga, f m*\ f 2 d\ gl [ d 1 ^ 

27r2 [mn) \dr^ ^ r dr) ^ 67 r^m* 

gu^ ^ (m* \ ( d^ 2 d\ g^g^ [ d ][ d ] 

■ S?‘Hm«) *3 + ;* ‘^■.('O + 5^ • (65) 




sea / \ e ^ f m* \ f 2 d\ . , 

Ppr,a(r) - -5^>‘i(^—j ( —+ - —j o(i- 




where 


dVl3\^) 

da 


dV)^\a) 

da 


i ~*~ ml ~*~ 2 ml I i ml ~*~ ml ) i ~*~ ml ~*~ 2 ml 


2 / &CT C(J^ ^ ^ 

\ ml ~*~ 2 ml I I ml ~*~ m 


b ca\ b'^ f. o •h\ ■hI 


-g <7 {Mm - ga(yf (^1 + 4 ln(l - + M^g^ - IM^gla 

+l?,MMgW - 


(9Z9)(c 


2 c{b + ca) 


{b + caf 


da 19277^ [(m^ + 6(T + Ica^) {ml + ba + ^ca'^y 
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In order to be able to calculate unclosed-shell nuclei, the occupation number have been 
explicitly indicated. Note that pQ®“(r) is a total derivative and thus the baryon number 
is conserved. The total energy of the system can be written as 


E — T'mft + Ah/, 


where 


£^mft = WaEo, - -J (fx -gatrps + -ba^ + g^cuoPo 

1 f 

+Y]^‘^oPo + 2^pRo,oPo,o + 4^^ ^opPofl + eAopprfi , 


AE = 


I d^x vYia) + vYia) + ^Z^^\a)iVaf 


^ln(^ 

47r2 \Mn 


67r2 Wiv 




The pairing energy and the center-of-mass correction to the total energy are taken into 
account as elucidated in Ref. [3]. The energy density of homogeneous nuclear matter 
can be obtained through reducing the above formulae. The compressibility at saturation 
density reads as 

9 m2 ° 3(/c^-|-m*2) 1/2 (p2, _|_ ^* 2 ^ 


X T(-) + + 


^PsT\ 


^=^0 (A;|,-|-m*2)i/2 m* 


where 


4 f^F 


(27r)3 7o (k2 + m*2)i/2- 

The double derivatives in the above expression can be easily computed. 


III. NUMERICAL RESULTS AND DISCUSSIONS 


Since the densities stemming from the Dirac sea are evaluated within the derivative 
expansion technique and expressed by means of the mean helds as well as their derivative 
terms, the wave functions of anti-nucleons are not involved when solving the meson-held 
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equations. The numerical procedure of the RHA is similar to the one currently used in 
the RMF [3], except that one more equation for the anti-nucleon is implemented. It goes 
as follows: in the n-th iteration step we have arrived at a set of wave functions of nucleons 
and mean helds. First, we calculate the densities contributed from the valence nucleons 
by adding up the nucleon wave functions, and the densities originated from the Dirac sea 
by evaluating the mean helds as well as their derivative terms. Second, we determine the 
meson helds by solving the Laplace equations of mesons. Third, we use the new meson 
helds to solve the Schrodinger-equivalent equation of the nucleon. Fourth, we compute 
the new single-particle energies of nucleons and determine the occupation numbers by 
adjusting a Fermi surface such that the particle number is conserved . This completes 
one iteration step. The iteration is continued until the binding energy is stable up to 
6 digits. Finally, we apply the known mean helds to solve the Schrodinger-equivalent 
equation of the anti-nucleon. The equation itself is solved iteratively. The obtained wave 
functions are used to calculate the single-particle energies of anti-nucleons. The space of 
anti-nucleons are truncated by the specihed principal and angular quantum numbers n 
and j with the guarantee that the calculated single-particle energies of anti-nucleons are 
converged when the truncated space is extended. We hnd that the results are insensitive 
to the exact values of n and j provided large enough numbers are given. We have used 
n = 4, j = 9 for n = 5, j = 11 for "‘^Ca; and n = 9, j = 19 for ^°®Pb. 

Table I 


The parameters of the model are determined in a least-square ht to the properties 
of spherical nuclei. The experimental values for the observables used in the £t are given 
in Table I. The second column gives the measured nuclear binding energies while the 
last two columns reflect the properties of nuclear shape. In model calculations one can 
extract the diffraction radius and surface thickness through analyzing the nuclear charge 
form factor, where the intrinsic nucleon form factors are buried in [3]. Here we just want 
to note that instead of the commonly used Sachs form factors [27] in the current htting 
we have applied recent parameterization of nucleon electromagnetic form factors based 


14 




on the Gari-Kriimpelmann model [28, 29, 30]. Specifically, we have taken the parameter 
set of GKex(02S) presented in Ref. [30]. 

In the fitting processes one can simply forget the anti-nucleon part since the vacuum 
contributions to the densities are calculated by means of the mean fields and their deriva¬ 
tive terms. Once the parameters are specified, we get a set of decided mean fields which 
are then applied to solve the eigenequation of the anti-nucleon. The eigenfunctions (the 
wave functions) and eigenvalues (the single-particle energies) of anti-nucleons are the fi¬ 
nal output of the model calculations. We intend, so to say, to predict the bound states 
of anti-nucleons through adjusting the model parameters to the bulk properties of finite 
nuclei. 

Gompared to the previous version of the RHA model now we have two more parameters 
foj and fp for the tensor couplings of vector mesons. The obtained parameters as well 
as the corresponding saturation properties are presented in Table II and denoted as the 
RHAT set. For the sake of comparison we give other two sets of parameters: the NLl set 
[3] of the RMF model under the no-sea approximation and the RHAl set [16] where the 
vacuum contributions have been taken into account but without tensor-coupling terms. 
One can see that after introducing the tensor couplings a large effective nucleon mass 
remains in the RHA model. The value of m*/M]^ ^ 0.8 is close to the one appeared in 
the Skyrme-force parameterizations for nonrelativistic approaches [19]. Although some 
rearrangements exist, generally speaking, the changes of parameters between the RHAT 
set and the RHAl set are not very significant except that the vector coupling strength g^ 
is enhanced. This causes a somewhat larger compressibility. 


Table H 


Table HI 


In Table HI the values of cire listed for the three cases and the deviations from 
the binging energy, diffraction radius and surface thickness are detailed. The spin-orbit 
splitting in and the shell fluctuation in ^‘^^Pb are provided together with the empirical 
values. Obviously, the NLl set performs a very good fit and describes the spin-orbit 
interaction satisfactorily. After including the vacuum effects the RHA model still reaches 
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a reasonable well fit to the bulk properties of spherical nuclei. The effect of tensor- 
coupling terms induces the spin-orbit force in the RHAT model to be one times larger 
than that in the RHAl model, and thus improves the total value. In addition, one 
can observe an interesting redistribution of deviations between the binding energy and 
surface thickness. In the RHAT set the x|; is suppressed substantially. It indicates an 
improved £t to the nuclear binding energies as clearly exhibited in Table IV. In the mean 
time, the contributions from the Dirac sea are enhanced evidently. 


Table IV 


Fig. 1 


The shell fluctuation can typically be expressed via the charge density in ^°®Pb. It has 
been known that both the relativistic and nonrelativistic mean-held theory overestimate 
5p by a factor of 3. From Table III one can find that the NLl set and the RHAT set share 
the same disease. However, the RHAl set reproduces the empirical value quite nicely. 
Further investigation is needed in order to clarify whether it comes out fortuitously. In 
Fig. 1 we depict the charge densities of three spherical nuclei reckoned with the RHAT set 
of parameters. The solid lines represent the experimental data. Except for the amplitude 
of the shell huctuations the overall results are in agreement with the data. 


Fig. 2 


Fig. 3 


The contributions of the vacuum to the scalar density and baryon density are shown 
in Fig. 2. The computations are performed with the RHAT set of parameters for ^°Ca. 
Noticeable inhuence from the Dirac sea can be found for the scalar density while the effect 
on the baryon density is negligible. We have solved the technical problem of huctuations 
on the met in the previous RHAl model through making spline extrapolation for the 
first several points of densities originated from the vacuum. Smooth curves for various 
densities in diherent nuclei considered in this work have been obtained. Fig. 3 shows the 
resultant scalar and vector potentials in for three models. Due to the vacuum ehects 
the potentials calculated with the RHA model are about half of that computed with the 
RMF model. After introducing the tensor-coupling terms, the RHAT set receives deeper 
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potentials compared to the RHAl set, reflecting the effect of parameter rearrangements. 
The enhancements are around 20 MeV for S and V in the center of the nucleus, which are 
nonnegligible on the scale of the nucleon central potential. Especially, the enhancements 
would be summed up for the anti-nucleon potential rather than cancel each other for the 
nucleon. 


Table V 


Table VI 


Fig. 4 


In Table V and Table VI we present the single-particle energies of protons (neutrons) 
and anti-protons (anti-neutrons) in three spherical nuclei of ^°Ca and ^°®Pb. The 
binding energies per nucleon and the rms charge radii are given too. The experimental 
data are taken from Ref. [31]. Both the relativistic mean-field theory (NLl) under the 
no-sea approximation and the relativistic Hartree approach (RHAl, RHAT) taking into 
account the vacuum contributions can reproduce the observed binding energies and rms 
charge radii quite well. With respect to the large error bars in measurements of the Is 
proton (neutron) levels, the results of all three sets of parameters coincide with the data. 
Because of the large effective nucleon mass, the spin-orbit splitting on the Ip levels is 
rather small in the RHAl model. The situation has been ameliorated conspicuously in 
the current RHAT model incorporating the tensor couplings for the oo- and p-meson. At 
the same time, a large m* stays unchanged. The experimental data for the anti-proton 
(anti-neutron) spectra in the vacuum are presently unavailable. The RMF model and 
the RHA model provide strikingly different predictions with a deviation of a factor of 
2, clearly demonstrating the importance of the Dirac sea effects. On the other hand, 
the anti-particle energies computed with the RHAT set of parameters are 20 - 30 MeV 
larger than that reckoned with the RHAl set, as can be anticipated from Fig. 3. The 
corresponding proton and anti-proton potentials in ^°®Pb are displayed in Fig. 4. The 
same features observed in the energy spectra are revealed once again. 

VI. SUMMARY AND OUTLOOK 

We have incorporated tensor couplings for the oo- and p-meson in a relativistic Hartree 
approach for finite nuclei. After refitting the parameters of the effective Lagrangian to the 
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bulk properties of spherical nuclei, the spin-orbit force has been enlarged by a factor of 2 
compared to the previous version of the RHA model without tensor-coupling terms, while 
a large effective nucleon mass remains. This improves the total value and brings the 
computed proton (neutron) spectra more closer to the data. The predicted anti-proton 
(anti-neutron) spectra in the vacuum are deepened about 20 - 30 MeV. One may argue 
that the vacuum may not be properly treated by the nucleon degrees of freedom, instead, 
a quark vacuum may be essential. However, if one speaks about observing a vacuum, what 
one actually means is to measure the response of the vacuum to the laboratory probes. 
In the environment of a hnite nucleus, nucleons are well established physical degrees of 
freedom, which should be the relevant degrees of freedom for describing the corresponding 
vacuum too. In the case that a QCD environment is involved, quark degrees of freedom 
may be necessary. 

In view that in the vacuum nucleons and anti-nucleons are always in the form of 
pairs, one of the promising ways to measure the anti-nucleon spectra in the vacuum of 
a nucleus is to knock out corresponding nucleons from the bound states emerging from 
the lower continuum. The incident particles can be photons, electrons or protons. The 
experimental searches for the vacuum structure have been discussed in detail in Refs. 
[2, 32]. The dynamical processes can be simulated by using the relativistic Boltzmann- 
Uehling-Uhlenbeck approach [33, 34, 35, 36] which is a microscopic transport model for 
single-particle distribution functions. Since the nucleons excited from the bound states 
of the Dirac sea have to overcome deep potentials in order to become real particles, 
the hnal particle spectra as well as the angular distributions should be different to that 
of nucleons originated from the Fermi sea. Exclusive analyses of observables from the 
hadron-nucleus reactions at the energy range of several GeV/c will exhibit the structure 
of quantum vacuum. Work on this aspect is in progress. If the energy of the incident 
particle is further increased, light nuclei could be directly excited from the vacuum when 
the correlation effect is taken into account. Relativistic quantum molecular dynamics 
model [37, 38, 39, 40] is a suitable starting point to study the relevant problems. 

It is straightforward to extend the present model to include the hyperon degrees of 
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freedom. Then one can apply it to investigate the properties of hypernuclei with the 
effects of quantum vacuum taken into account. As a hrst step we consider the single-A 
and double-A hypernuclei. Through systematically studying the A and nucleon spectra 
in the Fermi sea and the anti-A and anti-nucleon spectra in the Dirac sea, one can extract 
important information for the hyperon interaction, which is an active topic of modern 
nuclear physics [41, 42, 43, 44]. Here the anti-A spectra in the vacuum act as further 
constraints to the effective interactions in addition to the usual considered hypernuclei 
observables. The similar procedure can be performed for the S- and S-hypernuclei. The 
detailed understanding of the hyperon-hyperon and hyperon-nucleon interactions in dense 
medium is fundamental for the study of strange particle production and strange particle 
flow in relativistic heavy-ion collisions [45, 46, 47] as well as the composition and structure 
of neutron stars in astrophysics [48]. 
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Table 1: The experimental values for the observables included in the £t, the binding 
energy Eb, diffraction radius R and surface thickness a. In the last line we also give the 
adopted errors AO^ for the fit. 



Eb (MeV) 

R (fm) 

a (fm) 

160 

-127.6 

2.777 

0.839 

^oCa 

-342.1 

3.845 

0.978 

^sCa 

-416.0 

3.964 

0.881 

58Ni 

-506.5 

4.356 

0.911 

90zr 

-783.9 

5.040 

0.957 

ii^Sn 

-988.7 

5.537 

0.947 

i24Sn 

-1050.0 

5.640 

0.908 

208pb 

-1636.4 

6.806 

0.900 

AOJOn 

0.2% 

0.5% 

1.5% 
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Table 2: Parameters of the RMF and the RHA models as well as the corresponding 
saturation properties. Mm and rup are fixed during the £t. 



NLl 

RHAl 

RHAT 

Mm (MeV) 

938.000 

938.000 

938.000 

rria (MeV) 

492.250 

458.000 

450.000 

(MeV) 

795.359 

816.508 

814.592 

TUp (MeV) 

763.000 

763.000 

763.000 

ga 

10.1377 

7.1031 

7.0899 

gcu 

13.2846 

8.8496 

9.2215 

Sp 

9.9514 

10.2070 

11.0023 

b (fm“^) 

24.3448 

24.0870 

18.9782 

c 

-217.5876 

-15.9936 

- 27.6894 

J'uj/Mm (fm) 

0.0 

0.0 

2.0618 

Jp/Mm (fm) 

0.0 

0.0 

45.3318 

Po (fm"3) 

0.1518 

0.1524 

0.1493 

E/A (MeV) 

-16.43 

-16.98 

- 16.76 

m*/MM 

0.572 

0.788 

0.779 

K (MeV) 

212 

294 

311 

04 (MeV) 

43.6 

40.4 

44.0 
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Table 3: The standard values of the parameterizations given in Table II. The different 
sources of deviations from the binding energy, diffraction radius and surface thickness are 
separated. Other observables like the spin-orbit splitting of the Ip level in for both 
protons {Sep) and neutrons {Scn), and the shell fluctuation in ^°®Pb {Sp) are presented 
too. 



xl 

Xr 


xL 

Sep (MeV) 

SeN (MeV) 

Sp (fm 

NLl 

21.96 

11.78 

32.28 

66.02 

5.99 

6.06 

-0.0070 

RHAl 

516.48 

39.14 

256.69 

812.31 

1.99 

2.00 

-0.0030 

RHAT 

88.53 

24.50 

444.86 

557.88 

3.96 

4.43 

-0.0067 

Exp. 





5.98 

6.07 

-0.0023 


Table 4: Experimental and theoretical binding energies per nucleon within the RHA 
model. The vacuum corrections are indicated explicitly. 


RHAl RHAT 



Exp. 

Theory 

Dirac Sea 

Theory 

Dirac Sea 

160 

-7.98 

-8.00 

1.37 

-7.94 

1.68 

^°Ca 

-8.55 

-8.73 

1.43 

-8.62 

1.74 

"^Ca 

-8.67 

-8.51 

1.39 

-8.61 

1.74 

58Ni 

-8.73 

-8.44 

1.44 

-8.62 

1.85 

90 zr 

-8.71 

-8.74 

1.42 

-8.78 

1.76 

ii^Sn 

-8.52 

-8.61 

1.39 

-8.52 

1.68 


-8.47 

-8.50 

1.34 

-8.49 

1.68 

208pb 

-7.87 

-7.93 

1.30 

-7.88 

1.62 
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Table 5: The single-particle energies of both protons and anti-protons as well as the 
binding energies per nucleon and the rms charge radii in ^°Ca and ^°®Pb. 



NLl 

RHAl 

RHAT 

Exp. 






E/A (MeV) 

8.00 

8.00 

7.94 

7.98 

Eh (fm) 
PROTONS 

2.73 

2.66 

2.64 

2.74 

lsi /2 (MeV) 

36.18 

30.68 

31.63 

40±8 

lp 3/2 (MeV) 

17.31 

15.23 

16.18 

18.4 

lpi /2 (MeV) 
ANTI-PRO. 

11.32 

13.24 

12.22 

12.1 

lsi /2 (MeV) 

674.11 

299.42 

328.55 


lp 3/2 (MeV) 

604.70 

258.40 

283.44 


lpi /2 (MeV) 

605.77 

258.93 

285.87 


4uca 





E/A (MeV) 

8.58 

8.73 

8.62 

8.55 

rch (fm) 
PROTONS 

3.48 

3.42 

3.41 

3.45 

lsi /2 (MeV) 

46.86 

36.58 

37.01 

50±11 

lp 3/2 (MeV) 

30.15 

25.32 

25.95 


lpi /2 (MeV) 
ANTI-PRO. 

25.11 

24.03 

23.63 

34±6 

lsi /2 (MeV) 

796.09 

339.83 

367.90 


lp 3/2 (MeV) 

706.36 

309.24 

332.10 


lpi /2 (MeV) 

707.86 

309.52 

333.37 


2U8pb 





E/A (MeV) 

7.89 

7.93 

7.88 

7.87 

rch (fm) 
PROTONS 

5.52 

5.49 

5.46 

5.50 

lsi /2 (MeV) 

50.41 

40.80 

41.74 


lp 3/2 (MeV) 

44.45 

36.45 

37.38 


lpi /2 (MeV) 
ANTI-PRO. 

43.75 

36.21 

37.18 


lsi /2 (MeV) 

717.01 

354.18 

377.37 


lp 3/2 (MeV) 

705.20 

344.48 

366.95 


lpi /2 (MeV) 

705.28 

344.52 

367.24 
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Table 6: The single-particle energies of both neutrons and anti-neutrons. 



NLI 

RHAl 

RHAT 

Exp. 

TOq 





NEUTRONS 





lsi /2 (MeV) 

40.21 

34.71 

35.78 

45.7 

lp3/2 (MeV) 

21.07 

19.04 

20.18 

21.8 

lpi /2 (MeV) 
ANTI-NEU. 

15.01 

17.05 

15.75 

15.7 

lsi /2 (MeV) 

667.93 

293.23 

322.47 


lp 3/2 (MeV) 

598.74 

252.48 

277.94 


lpi /2 (MeV) 

599.74 

252.97 

279.22 







NEUTRONS 





lsi /2 (MeV) 

54.85 

44.48 

44.98 


lp3/2 (MeV) 

37.79 

32.98 

33.83 


lpi /2 (MeV) 
ANTI-NEU. 

32.73 

31.71 

30.99 


lsi /2 (MeV) 

783.87 

327.96 

355.70 


lp3/2 (MeV) 

694.80 

298.04 

321.07 


lpi /2 (MeV) 

696.18 

298.26 

322.15 


2U8pb 





NEUTRONS 





lsi /2 (MeV) 

58.97 

47.40 

46.70 


lp3/2 (MeV) 

52.44 

42.66 

42.31 


lpi /2 (MeV) 
ANTI-NEU. 

51.82 

42.45 

41.64 


Isi /2 (MeV) 

678.23 

313.18 

334.39 


lp3/2 (MeV) 

667.70 

304.61 

325.41 


lpi /2 (MeV) 

667.73 

304.61 

325.28 
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Figure 2: The scalar density and baryon density in ^‘^Ca. Dashed lines denote the con¬ 
tributions of valence nucleons, dotted lines represent the Dirac-sea effects and solid lines 
give the total resnlts. 
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Figure 3: The scalar potential from the a-meson exchange and the vector potential from 
the cn-meson exchange in Different curves are related to different sets of parameters 
as indicated in the hgure. 
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Figure 4: The potentials of the proton and the anti-proton in computed with 

different sets of parameters as indicated in the hgure. 
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